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Abstract 

Radiation emitted by unpolarized high-energy electrons penetrating crystals may be 
linearly polarized. This occurs when the particle velocity makes an angle, with respect 
to some major crystal axis, being sufficiently larger than the axial-channelling angle. 
For such orientation, a complete description of spectral and polarization characteristics 
of the radiation is derived. At planar channelling, a non-perturbative contribution to 
the probability of the process appears caused by the plane field, and we must solve 
exactly a one - dimensional mechanical problem. For that, the approximate form of 
the actual plane potential is suggested which provides a precise fit for any crystal plane 
and an analytical solution to the motion problem. In a practical case, we must consider 
electron-photon showers developing in sufficiently thick crystals. For the first time, this 
development is described taking into account the polarization of photons. We discuss 
qualitative features of the phenomenon, present results of numerical calculations for 
thin and thick crystals, and evaluate the possibility of the use of differently oriented 
crystals in a polarized hard photon source. 
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1 Introduction 



High-energy electrons emit the specific radiation while penetrating through single crystals. 
The shape of the spectrum and the intensity of this radiation depend on the electron energy 
and crystal orientation. This phenomenon has been widely investigated theoretically and 
experimentally ( see e.g., £Q and references cited there). Photon beams obtained in such 
a way may be used, for example, in photo-production experiments of different kind. In 
particular, they may serve as 7 - beams in 7 7-colliders. 

In this paper, we focus on the polarization properties of emitted photons as the polariza- 
tion offers additional opportunities in the experimental study. Initial electrons are assumed 
to be unpolarized. Then, using crystals, only linear polarization may be obtained. We start 
with the well known formula derived by means of the quasi-classical operator method ( see 
e.g., Remember, that within this method the probabilities of QED processes may be 
expressed by way of classical trajectories of charged particles, provided that the external field 
involved satisfies some conditions specified in [2j. Fortunately, the electric field of a crystal 
responsible for the coherent processes satisfies these conditions. In this stage, we have to 
solve a two - dimensional mechanical problem. 

Further consideration exploits the approach developed in As shown in [3], the me- 
chanical problem mentioned becomes essentially one - dimensional for the angles of incidence( 
w.r.t. some major crystal axis), i9o ; which satisfy the condition 6 sp < $0 <C 1. The magnitude 
of the angle sp depends on the electron energy and the axis chosen. It was accurately deter- 
mined in [3j being sufficiently larger than the characteristic angle for axial channelling . Note 
that the region $0 < 9 sp is not interesting here, as we do not expect any polarization of the 
radiation from unpolarized electrons at axial channelling. In the angular region of interest, a 
finite set of the most strong planes containing the axis is important. Moreover, the azimuthal 
angle domains, where the influence of the planes on a motion should be taken into account 
exactly, do not overlap. Then we have for the transverse (w.r.t. the axis ) velocity of an 
electron: v ~ v s i ow + Vf ast . In this expansion, the first term represents the exact solution 
for the transverse velocity in the field of the only one plane of the set. The second term 
gives the contribution of other planes which can be taken into account using a perturbation 
theory and, additionally, the rectilinear trajectory approximation. Note that periods of the 
one - dimensional motion ("slow") are noticeably larger then those in the perturbation term 
("fast"). The slow component is present when the velocity is almost aligned with one of the 
planes. It disappears when the velocity of a particle forms sufficiently large angles with each 
plane of the set. In this case, the perturbation theory in the whole crystal potential is appli- 
cable being the essence of the so called coherent bremsstrahlung (CB) theory. It is clear that 
the transverse motion of a particle is two - dimensional but the component which needs the 
exact calculation is one - dimensional. To find the latter, we propose the new approximate 
form of the actual plane potential which provides the precise fit for any crystal plane and 
allows one to find an analytic expression for the trajectory. That is done in Appendix A 
where also the explicit form of the velocity Fourier-transform is obtained. 

In Section 2, general expressions are derived giving the instantaneous (probabilities per 
unite time or length) characteristics of the radiation. Then they are analyzed and simplified. 
We discuss the qualitative features of the phenomenon and present results of numerical 
calculations for thin crystals. It turns out, that, for the two types of the azimuthal orientation 
mentioned above, not only spectra but also polarization distributions are utterly different. 
In a practical use, sufficiently thick crystals are needed to get a noticeable yield. In this 
case, we can not neglect the multiple photon emission, their absorbtion due to the e + e~- 



1 



pair production, and the radiative energy loss of charged particles. In other words, for thick 
crystals, we must consider the e + e~7-shower. Such a consideration is performed in Section 3, 
using the formulas obtained in Section 2 along with those describing the e + e~-pair production 
by polarized photons. So, we describe the shower development, taking into account the 
polarization of photons for both basic QED-processes involved. Some results of Monte-Carlo 
simulations are presented for settings used in NA43 ( see and NA59 ( see, e.g., 
experiments at CERN. The consideration performed allows one to decide for the optimal 
crystal type, thickness, and orientation providing necessary characteristics of a polarized 
photon source. 



2 Radiation at fixed particle energy 

Let us start with the well known formula ( see, e.g., Eq.(16.7) in pQ ) which describes 
spectral, angular, and polarization distributions of photons emitted by unpolarized electrons 
(positrons) at given motion: 

a d^k s f s 

dw^ = — — — / dt 1 dt 2 L(t 1 ,t 2 ) expU—^k^x - x 2 )} , (1) 

{Ajc) z to £ J £ 

where a = 1/137 , k^ = {u),k) is the photon momentum, x% 2 = 2^(^1,2), x^(t) = (t,x(t) , 
e' = e — uj, e is the electron energy, and 

L(t ls f 2 ) = (e* Vl )(ev 2 )(if(£) + 2) + [(e*e)(«i« 2 - 1 + 7 ~ 2 ) - (e*v 2 )(e Vl )](<p(£) - 2). (2) 

Here V\ )2 = ^(^1,2) is the electron velocity on the classical trajectory x(t) , 7 = e/m, m is the 
electron mass, ip(e) = e/e' + e'/e. If we set ip(e) = 2 {uj = 0)in the expression for L{t\,t 2 ), 
the formula will describe the radiation from the scalar (zero-spin) particle ( see e.g., part 
9 in |H] or P ). This substitution rule holds for all subsequent formulas in this Section. 

Remember that the vectors e in Eq.(J2J) correspond to the polarization which would be 
measured by some detector and, thereby, have nothing to do with the polarization of emitted 
photons. The latter is described by the matrix, dwu in the contraction dw 1 = eie* k dwki- 
We introduce the vector, £ (|£| = 1), which describes the analyzing ability of a detector, by 
eie* k = (1 + £er)j fc /2 , where cr are the Pauli matrices. The matrix dwui can be presented in 
the same form: dwki = (A + Ba) ki /2 . Then we have 

1 A B 

dty 7 = -{A + B$,) = — (1 + ry£) , rj = — , A = dw u + dw 22 , 

(3) 

Bi = dw 12 + dw 2 i , B 2 = i(dw± 2 - dw 21 ) , B 3 = dw u - dw 22 . 

In this equation, the Stokes vector of the radiation, rj, is defined. The quantity A gives the 
value of dwry summed up over polarizations. 

We assume further that the angular divergence of the electron beam is small enough. 
Then the same is true for arising photon beam provided that electron energy is sufficiently 
large (e ^> m), as the emission angle of a photon (w.r.t. the particle velocity ) is typically of 
the order of 7 _1 <C 1. In this case, we can choose some axis (z-axis) such that the momenta 
of all charged particles and photons make small angles,??^, w.r.t. this axis: flax *C 1. Let a± 
denotes the component of an arbitrary vector a transverse to the axis chosen. Using relations 

v z ~ 1 - ^{T 2 + «1) , e 2 = -(k ± e ± )/k z ~ (n±e±) , n z ~ 1 - ^n 2 ± , n = k/u , 
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we can remove z-components of all the vectors from Eq. (|TJ) . As a result, we have 
aujduod 2 n±_ e 



(4tt) 2 : 



t 2 



-J dt\dt 2 L{t\, t 2 ) exp(— iD) 



D = %pj dt[T 2 + (n±-v ± (t)f 



L{ti,t 2 ) = (e* L ,v 1 _ L -n_ L )(e_ L ,v 2 ±-n_ L )(<f(e) +2) 



(4) 



- [(el, w 2 ± - n±)(e±,v 11 - n±) + -{e^e^v^ - v 1± ) 2 ](ip(e) - 2) . 

No restrictive conditions (like (ne) = 0) are imposed on the vectors e±. This allows one, 
first, to utilize any fixed basis for the description of the photon polarization regardless of 
n. In what follows, we use the Cartesian basis (e x , e y , e z ), where e z is directed along the 
nearest (making a small angle w.r.t. the particle velocity ) major crystal axis , e x is within 
some crystal plane containing this axis, and e y is perpendicular to this plane. Secondly, the 
vectors e± can be considered constant when integrating in (J3J) over photon emission angle ( 
over d 2 n± ). Performing this integration, we obtain the matrix, w^, describing the spectral 
distribution of a radiation: 



Lij 



iaduj 

47T'7 2 



dt 



dr 



T-iO 



exp[— i\r( 1 + p{t, r))] L 



(5) 
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\^)-2)( 92 - 9l f + w ^- Q)i 



+ 



9%i9\i ~ 9ij9\i - \92i9\j + 92j9u 



where t = (ti +t 2 )/2 , r = t 2 -ti, A = um 2 / (2ee f ) , g^ 2 = g(ti,2) , and 

t/2 t/2 



9{s 
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v ± is 



1 

T 



p(V 



1 

T 



dx g 2 (t + x). 



dx v±(t + x) 

-t/2 -t/2 

In further consideration, we use the approach developed in [3]. Recollect that the trans- 
verse velocity can be represented as the sum: ~ vf° w + v^ st ', where v^ st given by 
Eq.(3) in [3] is characterized by small amplitudes and large frequencies. On the contrary, 
large amplitudes and small frequencies are inherent in the term vf° w which corresponds 
to the one - dimensional motion in the field of a plane. Its explicit form reads vf° w (t) = 
e y S y^ j v n exp(inu Q t) . Here u Q is the frequency of this motion, and v n is the velocity Fourier 

n 

transform calculated in Appendix A. Correspondingly, the quantity g(s) in (JHJ) turns into 
the sum: g(s) = l(s) + w(s) where 

sin(nu;oT/2) 



l(s) 



w(s) 



itiujqs 



irudot 



mq\\ 



tiwqt/2 



sin(g||T/2) 
9IF/2 



(6) 



J(q±po) 
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Here q are discrete reciprocal lattice vectors, G(q) are coefficients in the fourier series pre- 
senting the crystal potential ( see Chapter 9 of [1^ for the explicit form of q and G(q) ), 
and gii = (q±v) . The tilde in the expression for w(s) means that the sum does not contain 
q± || e y which just form the plane potential. 

In principle, using (J5jl. one can perform the integration over t and r straight in ( | 5 jl . 
However, as we have already obtained the quantity w(s) by means of some perturbation 
procedure, it is more consequently to continue in the same way. So, as in [3], we expand the 
exponential function in (jSj) in w(s) , keeping quadratic terms, and obtain 



dw-f - - — - 

47T7 2 



oo oo g 

iadu f f dr , . , ... v-^ r n ) 

dt I exp[-t\r( 1 + p s i ow {t, t) )J > Q, 

/ T — ZU * — * 

^ n=l 

-oo — oo 



(7) 



where p s iow(t, r) = r 1 J^/ 2 cfo l 2 (t + x) and matrices are 



C 



(i) 



A(r-iO) 4 



+ 7 (^(e)-2)(« 2 -«i 



(2) 



E 



l-/ 2 (C) + (^)-2)sin 2 C 



/( 2 +) (C) + /( 2 -)(C) 



c 



(3) 



2 ^E 



J, 



J, 



Kl 



T 



no 



-c. 



(1) 



A(r - iO) u J 



+ zA^fC^ + $<{ z<^(e) - 2)(i 2 - I ls i/) sinC 



Vie) 



+ (iii^- + hjVi)f (-)(C) + (fe^ + %^)/(+)(C) + (iii^i - hjVi)! (-)(() 



{kiVj -hjVi)f (+){() 



t/2 
-r/2 



/(C) 



/(±) ( x ) =e ± --/(x), »/ 



9± 



C = 9||T/2. 



Using Eqs. 0,©, we obtain the expression for the probability of a photon detection irre- 
spective of its polarization, dw^ np = dw^S^ , wich coincides with Eq.(4) in [3]. 

Let us dwell on the plane field contribution (PFC) to the radiation. This is given by 



the term in J7|) which is proportional to C^- . To make estimates, we should recollect some 
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properties of the quantity p s iow(t, t) which enters into the phase, ^(r) = Ar( 1 + p s iow(t, t)). 
So, p s iow(t,r) is the even, positive, and monotonically increasing function of r . It satisfies 
the inequality p s iow(t,T) < p(oo) < p c = U p ie/m 2 where U p i is the potential well depth of 
the plane and p(oo) is the limiting value of p s iow(t, T ) as |r| — > oo . Using (JHJ), we can easily 
estimate the behavior of p s i ow (t, r) at large and small r 



< P?, > - < Py > 2 



p(u; r < 1) 
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(iv y (t)r) 2 , 



(9) 
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where < . . . > means time averaging, p y and v y are correspondingly y-components of the 
momentum and acceleration of a charged particle. Let $0 and 0o be the polar and azimuthal 
angles of incidence counted off from the axis and plane correspondingly. Then the velocity 
of a particle makes the angle, ip, w.r.t. the plane, at that, simp = sin^ o sin0o or, as $0 *C 1 
, we have ip ~ ^osin0 o . At channelling ( ip < 9 p i = {2U p i/ e) 1 ^ 2 ) , we have p(oo) < p c . 
At above-barrier motion ( ip > 8 p i), the quantity p(oo) decreases fast being of the order of 
0.l(U p i/miP) 2 for if) > pl (see, e.g., p.430 in 0). 

For p(oo) C 1 , we can expand the exponential function in (J2J in powers of p s iow(t, T ) 
retaining only the linear term ( dipole approximation). In this case, the slow-fast interference 
term, , gives higher-order corrections and should be neglected, while after averaging 
over time t reads 



C\)\p « 1) = J2 l7^| 2 {^[l-/ 2 (C) + (^(^)-2) sin 2 c] -e yi e yj [f 2 +) (() +/f_ ) (C)] } 



(10) 



where £ = uujqt/2 , = 5 2 i- This expression has the same structure as ( cf. (JBJ)) since 
both were obtained by means of the perturbation theory. For p c 1 , Eq. lfTUj) is valid at 
any value of the angle ip. In the opposite case when p c > 1 , it holds for ip ^> 9 p i. The latter 
condition provides also the applicability of the rectilinear trajectory approximation in the 
calculation of v n and uj . So, as we have checked using the explicit form of v n and Uq ( see 
Appendix A ), for ip ^> 9 p i , and are essentially the same except that the summation 

is performed in over the subset of q± parallel to e y , while in Cj-j the complementary 
subset is used. As a result, when the particle velocity is well off all major planes, the spectrum 
and polarization of a radiation is described by the term C^ 2 if we extend the summation in it 
over all q±. Then it reproduces results of the so-called coherent bremsstrahlung (CB) theory. 
In [7]( see Eqs. (4. 8), (1.5), the correlation was explained of this theory with the Compton 
scattering (off electrons) of the equivalent photons representing in a proper reference frame 
the periodic crystal field. 

For p(oo) ^> 1 , the constant field approximation (CFA) is widely used. It can be 
applied when the magnitude of the external ( not uniform ) field is almost constant on the 
particle trajectory during the radiation formation time, Tf . The latter is determined by 
the condition ^(t/) < 1 . In our case, the applicability condition of CFA reads u Tf <^ 1 . 
As Tf depends on photon energy, uj , this condition can not be fulfilled everywhere in the 
radiation spectrum. In particular, it is violated for uj < 27 2 o; /p(oo). The power spectrum 
is maximum at u = u)/(e — u) ~ x — l 2 v{t)/m . Using Eq. Q, we find at u > x f° r the 
formation time Tf ~ 2u /('-fv(t) . Numerical estimates of the latter expression show that, 
even in this part of the spectrum, the CFA becomes valid at rather high electron energy. For 
example, it happens at e ^> 100 GeV for (llO)-plane of diamond and silicon, and at e 3> 10 
GeV for the same plane of tungsten. Using CFA, the probability, dw^, was calculated in [E]. 
We emphasize that no approximations are used in the present paper at the PFC calculation. 

However, the results of [H] are useful here since they allow us to check the procedure of the 

(3) 

exact calculation and estimate the magnitude of the interference term, , as well as the 

(2) 

influence of the slow motion on the CB term . This influence is due to the presence 
of p s iow(t, T ) in the phase ^(r) . According to [H|, the scale of the effect is given by the 
parameter p = xl s (.1\\) where s(q») = 2e\q\\\/m 2 is a conventional kinematic parameter in 
the Compton scattering . Remember, that the power spectrum of the latter is maximum just 
at the kinematic boundary, u = s(qn). So that, the plane field and CB contributions to the 
radiation intensity are well separated for p 1 . The parameter p is essentially the deviation 
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of the particle velocity due to the external field action during the hard-photon formation time, 
T h ~ kill -1 j measured in units of the characteristic emission angle 9 p h = 7 _1 - For given 
plane and substance, this parameter depends only on the angle of incidence, $ > namely, 
\i = C^/$o(mi&d). Using the Tables 9.2 and 15.1 in jl] , we obtain for < 110 > planes, 
(diamond) ~ C M (silicon)~ 0.01 , and (^(tungsten) ~ 0.06 . For weaker planes of the same 
crystal, C M is smaller being proportional to the magnitude of the plane electric field. From 
Fig.l in jSj, the influence of the field on the Compton scattering can be neglected if \x does 
not exceed several hundredth. In turn, the interference term provides corrections to PFC of 
the order of pf ast ~ ^2\G(q±)q±/ (mq\\) | 2 which, by definition , is much smaller than unity. 
Assuming /i < 0.05 , we obtain from a simple but still rather accurate expression for the 
photon emission probability, dW, per unite length (time) as well as the spectral distribution 
of polarization. According to the analysis performed in jH], the accuracy of 5 % or better is 
expected. Let us present dW in the form introduced in as just this form will be used in 
further calculations: 



dW 

did 



+ 



did 



duo 



did 



where the superscripts C, F are correspondingly for the CB contribution and PFC. For the 
latter (plane field contribution), we have 



OO 7T 



a Jdx^ F \b^) ] b[ F) =bi F) = 0,bi F) = (Dj-Dl)sm* 

cos(s5) — cos \I> 



o o 



(cp(e) - l)D( - D{ 



sin* + 



s5 



* = s5(l + D 3 ) 



:i2) 



D\ = 27 v n sin(ns) sin(nx) , D 2 = 27 v n g{ns) cos(nx) , 

n=l n=l 

OO 

27 2 v n v m {[f((n + m)s) — f(ns)f(ms)} cos((n + m)x) + (m — > — m)} 
/( X ) = _E^ > g(x) = cosx - f{x) , 



D 



n,m=l 



m 2 u 



euj 



where the integration over x = id^t corresponds to time averaging, the variable s = ojqt/2 
is used instead of r . When the influence of the plane field on the Compton scattering is 
neglected ( the factor exp(— iXrp^owiti r )) is omitted in the term oc in Q), the integrals 
over r are easily taken and we obtain for the CB contribution: 



A (C) B {C) 
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E 



G(q±)q± 



mq\\ 



u 



ijjm 



M\\) 2ee'|g||| 



(13) 



,(C) 



W u(o)_ b (c) 



1 "3 



where 9(1 — j3) is the step function: 9(x) = 1 for x > and 9{x) = for x < , v is defined 
in (H. 
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Figure 1: (a): intensities xd 2 N 1 /dxdl at e = 180 GeV for (llO)-plane of silicon at z = —0.25 
(solid), at z — 1 (dash-dotted); the same in diamond at e = 180 GeV , z = —0.25 for (1 10)- 
plane (dashed) and for (OOl)-plane (dotted), (b): third component of the Stokes vector for 
these conditions except z = 1 ; z = e±/U p i , x — uj/e. 

Using formulas obtained, we present now some examples illustrating the characteristics of 
a radiation. Let us start with the PFC described by (|12|). In FigQ, the radiation intensities 
(a) and polarizations (b) are plotted as functions of x = uj/e. Remember, that the quantities 
v n and ujo in (fT2"|) depend on the integral of motion, e± and so does the radiation. From 
Fig. , the radiation at channelling (z = e±/U p i < 0) is softer and more intensive than 
that at above-barrier motion {z > 0). At given z ,the intensity is smaller for a weaker plane 
(cf. curves for diamond). As expected, the polarization at such one - dimensional motion is 
directed perpendicular to the plane (r/ 3 < , r]\ = 0). The polarization degree is rather high 
and does not reveal a sharp dependence on x. 

In FigJ21 , the radiation intensity, probability, and polarization are shown for the angles 
of incidence i?o = 5 mrad and 0o — 0.036 when ip = 180/irad. For the (llO)-plane of Si at 
e = 180 GeV, the channelling angle is of Q v \ ~ 15/xrad,i-e., ip 3> 9 p i (z ~ 140). Then, as 
explained above, the first term (PFC) in (jllj) has the same form as the second one and we are 
dealing with pure CB. This term is dominant when, as in our example, (f>o -C 1. That is due 
to the relative smallness of gP from the first subset (roughly gP ~ gP (f> ). As a result, the 
second term has a much smaller amplitude ( g7 2 in a partial flux of equivalent photons ) and 
much higher photon frequencies (/3 oc q^ 1 in (fT3j)). So, it can be neglected for (p <^ 1 and 
the expression (JTHJ) is reduced to a one - dimensional sum. Within this accuracy, we obtain 
for r/ 3 at the maximum of the first harmonic, i.e., at u = s : 

•(•-^--{Ei^rMEi^r^^-io-s)]}. <"» 

n=l n=l ' 

where s = 2eqip/m 2 , q = 2n/d p i , and d v \ is the inter-planar distance. The contribution of 
the first harmonic (n = 1) to r/ 3 at u = s is rj^ 

L > = -[1 + s 2 /(l + s)/2) l . Worthy to note 
that 773^ is independent of G(q). In our case, when s 2 /(l + s)/2 ~ 0.49 , it overestimates 
the exact value (|T4*jl by 8%. Recollect now, that according to E Sj each equivalent photon 
is completely linearly polarized along its q±. In the first subset, all such photons have q± 
perpendicular to the plane. Thus, the whole equivalent photon beam produced by this subset 
is completely linearly polarized leading to the polarization of emitted radiation perpendicular 
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to the plane (773 < , T]i = 0). Actually, the second term in dW/<ko (see (JUJ) has not been 
neglected in our calculation giving at this orientation non - vanishing but extremely small 
value of r]\. To give a glimpse of a role of the electron spin, the intensity and polarization 




Figure 2: (a): solid curves present intensity xd 2 N^/dxdl and probability d 2 N 1 / dxdl (upper 
curve) at e — 180 GeV of CB in silicon ( 5 mrad off the < 001 >-axis and 180 /irad off the 
(llO)-plane ), dotted curve presents intensity for a "scalar" electron ; (b): third component 
of the Stokes vector at the same settings 

from a " scalar" (zero spin) electron are presented in Fig|2^ dotted curves) . The spin terms in 
(J2J) and subsequent formulas are proportional to ip{e) — 2 = u 2 / {se'). Therefore a difference 
in radiation characteristics becomes observable at sufficiently large photon energies ( for 
x > 0.3 in Fig|5J). The peaks do not move being determined, at given orientation, solely by 
the particle energy, e . Due to the absence of the spin terms, the radiation from a "scalar" 
electron is less intensive and, because of that, has a higher polarization. Really, the quantity 
B^ in ()13j) is independent of the particle spin, while A^ c \ which is proportional to the 
intensity, appears in the denominator of the equation defining a polarization ( 77 = B/A). 
In particular, the component 773 (u = s) for a "scalar" electron is given by Eq. (jl4j) if we omit 
the item s 2 /(l + s)/2 in the denominator. 

When both terms in (JTTJl contribute to the radiation, a corresponding alignment is some- 
times called string - of - string (SOS) orientation, since at such an alignment particles traverse 
axes (strings) forming the plane. In this case, the emission of hard photons is described by 
the second term in (jllj) being CB by nature. A difference of spectra and polarizations for 
photons emitted by means of the same physical mechanism (CB) is completely due to that 
of equivalent photon fluxes at different orientations. The quantity which determines the 
shape of a spectrum is independent of q 2 at SOS: q\\ = qx'&Q. Then the summation over qi 
in — Yl ( see corresponds to the splitting of the total flux into subsets of equiva- 

lent photons having the same q». Remember, that for x < x max (x max marks the first peak 
position ) the radiation is described within a high accuracy by only one subset having the 
minimal value of \q»\. This value is provided by |<7i| m j n = 2n/d ax where d ax is the distance 

between axes forming the plane. As the quantity 63 in ()13|) is proportional to q\ — g|, only 
first (q 2 = 0) term in the sum over q 2 is positive for the main (l^l = |5i| m » n ) subset. When q 2 
increases, the magnitude of negative terms diminishes. Their sum, however, cancels almost 
perfectly the first term. This results in a small magnitude of 773 seen in FigOl For the next 
(|gi| = 2\qi\ min ) subset, already two first terms are positive which leads to a positive total 
sum over q 2 . As a result, the polarization is somewhat higher for x > x max being parallel to 
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Figure 3: (a): Contribution of the CB-like term ( Eq. fjlljl at SOS orientation ( 0.3 mrad off 
the < 001 >-axis in the (llO)-plane ). Intensities and probabilities in diamond at e — 150 
GeV (solid ), and in silicon at e = 180 GeV (dotted) ;(b): third component of the Stokes 
vector at the same settings. 



the plane (773 > 0). Qualitatively, the equivalent photon beam at SOS orientation turns out 
to be almost unpolarized in contrast to the above example of pure CB. The component 771 
vanishes at SOS orientation. This can be easily verified if we, e.g., change the sign of q 2 in 

the sum ^2 in (fT3*j) . 

92 

The total ( integrated over u ) probabilities, W* ot , are typically the order of magnitude 
larger at SOS orientation mainly due to the PFC. More precisely, we have in the above 
examples for silicon W* ot (CB) ~ 1.2cm -1 , while at SOS orientation the PFC and CB-like 
terms give correspondingly 42.9cm _1 and 2.2cm _1 . 

To describe a shower development, the probability of e + e~-pair production by a photon, 
dW e , is needed as well. Using Eqs.(3.12),(3.25) in Ij, we obtain, first, the expression analo- 
gous to (jlj where the summation over positron final states and electron spin states has been 
performed 

a d^p s f s 

- / dt 1 dt 2 L e (ti,t 2 )exp[-i—(k,x 1 -x 2 )]. (15) 



dw P 



(2?r) 2 uo 2e' 

Here p and e are the electron momentum and energy, e' = uj — e, and 

L e (t u t 2 ) = {e*v 1 ){ev 2 ){ip{s) - 2) + [(e*e)(«i« 2 - 1 + 7 ~ 2 ) - (e*r 2 )(e«i)](^(e) + 2). (16) 

Note that the quantities L e (ti,t 2 ) in (|T6*|) and L(t\,t 2 ) in (J2J) turn into each other if we change 
(<p(e) — 2) < — > (ip(s) + 2). Further consideration of the pair production may be performed 
using the same approach and approximations as those applied above to the photon emission 
problem. Here we give explicitly only a perturbation ( CB-like ) term in the expression for 
dW e in the form of (fTTjl 

dW e _ dwP dwP dW^ F) _ iCjF) {G)F) 

where 77 describes the photon polarization, and 

(^,B f »)^ E |^| 2 ( a f», 6 f.), (1 -« ; P= ^ y 
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(18) 

aP = \<p(e) + 13(1 - (3) , b$> = -0*w , b£> = , bg> = -\^\ ~ 4) ■ 

Note that this result may be obtained from (|13|) by means of the substitution mentioned 
above and evident change in the common multiplier. The plane field contribution to the 
pair production probability ( first term in fTTj) ) was investigated in 9\ where the CFA was 
used. Though the applicability of this approximation is questionable ( see discussion above 
), we used the results of [9 J as a rough estimate and found that PFC to the pair production 
probability should be neglected under our conditions. 

3 Radiation from thick crystals 

As long as the crystal thickness, L, satisfies the condition iV 7 ~ W^ ot L <C 1 (thin crystal), 
the radiation emitted is described by formulas obtained in the previous Section. For such 
thicknesses, the relative energy loss, Ae/eo, is even smaller than the number, iV 7 , of photons 
emitted. Since the total probability, W" 7 0< , depends on the initial electron energy e and 
crystal orientation, the same sample may prove to be thin or thick (iV 7 ~ 1) depending on 
settings. At the noticeable (iV 7 > 1) yield, an alteration of the particle energy can no longer 
be neglected, several photons are emitted and the electron-photon shower develops. 

The main processes taken into account in our simulation of the e + e~7-shower development 
are: i) emission of photons due to the coherent and incoherent mechanisms, ii) absorbtion 
of photons due to the e + e~-pair production by both mechanisms, iii) multiple scattering 
of electrons and positrons. Many-dimensional maps of probabilities were created describing 
the photon emission and pair production depending on the energy, momentum direction, and 
polarization. In other words, thousands distributions like those shown in Figs llJ2J3*l have been 
obtained while a calculation of the PFC was the most arduous task. Different mechanisms 
were simulated as independent ones. In particular, each photon emitted by the coherent 
mechanism was provided with the polarization according to Eq.()ll|) and was unpolarized 
when emitted incoherently. The angular divergence of the initial electron beam was also 
taken into account. The values for this divergence of 30 /xrad and 50 /xrad used in our 
calculations correspond to the experimental conditions of [Sj and [3] as do also the initial 
energies and angles of incidence. Note that the same settings were used in above examples 
(see Figs|2J HU) illustrating instantaneous characteristics of a radiation. So, we can compare 
the outputs from thin and thick crystals. 

One must distinguish the true power spectrum of a radiation from that of energy losses. 
The latter is observed when a detector (e.g., a calorimeter) sums up over the energies of all 
photons emitted by one electron. These spectra coincide in the limit of vanishing crystal 
thickness when the multiplicity (a number of photons at given nonzero total energy loss) 
tends to unity. In our case, they are very different, especially for SOS orientation (cf. curves 
1 and 2 in FigEk) when the PFC is dominant leading to rather soft power spectrum. Note 
that due to the same mechanism (PFC) of photon emission, the multiplicity is very high just 
for SOS orientation (cf. curves 1 and 3 in Fig. U>). Remember that the multiplicity increases 
when the lower energy boundary of photon recorded , u> t h, decreases. Results presented in 
FigEb were obtained at u th = 1 GeV. In applications, power spectra are less interesting than 
the distributions in number of photons (spectra). Such spectra are presented in FigOJj in the 
form of effective probabilities capable of direct comparison with probabilities shown in FigsEl 
EJ The shower spectra are significantly softer than initial ones due to the decrease of the 
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Figure 4: Yield from 1.5 cm-thick silicon crystal at (center of beam) settings as in Figs J2J3| (a): 
Power and energy loss spectra for SOS (1,2) and CB (3,4); (b): spectra (effective probabilities) 
for SOS (dotted) and CB (solid). 

mean energy of charged particles with the increasing depth, and to the incoherent mechanism 
action. Recollect, that the parameter |cj||| which determines the position of hard peaks in the 
instantaneous spectrum depends on the current energy and velocity direction. The latter also 
changes in thick crystals mainly due to the multiple scattering. Note that for conditions of 
FigHJ the mean-square scattering angle ( at the initial energy ) is about 50 /irad being larger 
than the angular divergence of the initial electron beam. As a result, a smearing of peaks 
takes place and sharp structures are not seen in shower spectra. All the factors mentioned 
affect the shape of polarization distributions shown in FigEK For CB ( see curve 3 in FigEK 




co (GeV) energy loss (GeV) 



Figure 5: (a): Polarization (7/3) at settings of FigjUfor SOS (1), CB (3), and for SOS from 
0.05 cm-thick diamond crystal (2) at settings of Fig.l in 4J( Eq = 150 GeV, i9n = 0.3 mrad 
|^| < 10/xrad) ; (b): multiplicities at these conditions for u t h = 1 GeV. 

), this distribution is maximal at u ~ 100 GeV. The shift left (by ~ 11 GeV as compared 
to FigEb) is due to the energy and angular spread mentioned above, while a diminution of 
the magnitude is mainly caused by the incoherent (unpolarized photons) contribution. For 
evident reasons, such changes are marked feebly for relatively thin diamond ( see curve 2 in 
FigEK ) crystal. 

Let us define an enhancement as bin-by-bin ratio of the yields from oriented to disoriented 
crystals. In the latter case, only incoherent (amorphous-like) mechanism is acting, when 
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dNj /dx ~ Q(x)/x with Q(x) being rather smooth function of x. Therefore, the ratio of 
numbers of photons ( probability enhancement ) is very like in the form to a power spectrum 
while the energy loss enhancement bears a strong resemblance to the energy loss spectrum 
shape (cf. corresponding curves in FigsJ3^ and )• Note that an enhancement increases 
with decreasing thickness L, other things being equal. This explains, along with a larger bare 
(at L — > 0) probability enhancement for diamond, the order of magnitude difference in these 
quantities for silicon and diamond crystals (cf. curves 1 in Figs. andEb )• Already from 




40 80 120 160 20 40 60 80 100 120 140 
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Figure 6: At settings of FigE3 (a): probability and energy loss enhancements in silicon for 
SOS (1,2) and CB (3,4); (b): the same for SOS in diamond 

Figs. 141 51 the SOS orientation looks less favorable than CB for the hard photon production. 
As explained above, this is due to the PFC which in itself is characterized by relatively 
soft spectra with large intensities and total probabilities. So that, the CB-like contribution 
providing a hard photon emission is suppressed at SOS, in particular, due to the energy loss 
via competing mechanisms. However, the thickness in above examples was chosen to optimize 




Figure 7: Number of photons per bin Ax = 0.02 in a diamond crystal at Eq = 250GeV for 
CB ( 5 mrad off the < 001 >-axis and 160 /xrad off the (llO)-plane ) at L — 0.1cm (1) and 
L = 0.5 cm (2); for SOS ( 325 /iiad off the < 110 >-axis in the (OOl)-plane ) at L = 0.04 cm 
(3) and L = 0.2 cm (4); x — tu/e . 
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the yield of CB at u ~ 100 GeV being not optimal for SOS. Additionally, the positions of 
peaks (cf. Figs. |21and|Hl) in the initial spectra were different for two orientations. Let us now 
compare the yield from a diamond crystal at Eq = 250 GeV for three different orientations 
characterized by the same peak position in instantaneous spectra, x max = 0.75. Those are 
two SOS orientations ( 325 /irad off the < 110 >-axis in the (OOl)-plane and 180 /irad off the 
< 001 >-axis in the (llO)-plane ) and CB ( 5 mrad off the < 001 >-axis and 160 /irad off the 
(llO)-plane ). From Fig. Q the weaker (OOl)-plane is preferable to the stronger (llO)-plane 
for the purpose of hard photon production. Additionally, stronger < 110 >-axes are involved 
in the first of the two SOS orientations increasing the CB-like contribution. As a result, 
the hard (x > 0.5) photon yield turned out to be higher for this orientation at any crystal 
thickness. For example, this is 1.5 times as large in the region x ~ 0.6 -j- 0.75 at L — 0.2 
cm. Just this orientation is compared to pure CB in Fig. [7| For CB, a local maximum in 
the shower spectra situated at x = 0.69 is seen even at larger L than those presented in 
Fig. [7| At SOS, such a maximum (at x = 0.71 ) is feebly marked already at L = 0.04 cm, 
so that the spectra are monotonically decreasing for larger L. However, an increase of the 
yield at x = 0.71 continues up to L = 0.3 cm where it amounts only to 3% of the yield at 
L = 0.2 cm. Thus, a saturation occurs in the hard part of the SOS-spectrum at thicknesses 
L ~ 0.2 cm. In our examples of CB, the thicknesses are chosen so as being noticeably smaller 
than those leading to the saturation of CB-spectra, they provide almost the same amount of 
hard photons as at SOS. Even under these conditions ( the hard photon yield of CB may be 
further increased ), CB-spectra would get the better in applications being on the whole much 
harder than SOS-spectra (cf. curves 1,3 and 2,4 in FigJZJ), and having a polarization ( about 
40% near x max ). Numerous relatively soft photons may provide rather severe background 
conditions, producing, e.g., e + e~-pairs or hadrons directly in a radiator. 

In conclusion, we hope that explicit formulas presented along with a qualitative analysis 
performed allow anyone to make own estimations of radiation characteristics at any orien- 
tation where polarized photons may be obtained from unpolarized electrons or positrons 
penetrating through single crystals. 



Appendix A 

Let x be the coordinate perpendicular to some system of crystal planes with the inter-planar 
distance, d p i. Then the periodic plane potential for electrons reads 

U e {x) = - V" G(nq) exp(iimy) , g=y^, V = ~r ■ 

dpi dpi 

n=— oo r r 

As the direct use of this potential is impossible in analytical calculations, several approximate 
forms ( see Chapter 9,15 in ) were suggested. Being not satisfied with previously used 
forms, we propose here the new one which provides a precise fit for any crystal plane potential 
and very simple expressions for the velocity Fourier-transforms. For < y < 1, it reads 

U e (y) = -f/ P 4^(y 1 -y)(l-a 1 y 2 )+^(y-y 1 )^(y 2 -y)(a 2 y 2 +6 2/ + C )+^(y- 2/2 )a3(l-y) 2 ] , (A.l) 

where 9(x) was defined in (|T3*j) being the step function, U p i is the potential well depth. The 
origin is set to the point where the potential is minimal {U e (0) = —U p i),i.e., just at the 
plane. Beyond the segment y G [0, 1], the values of U e (y) may be obtained from (|A.1|) using 
evident symmetry and periodicity conditions. For positrons, we also choose the origin at the 
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point where the corresponding potential, U p (y), is minimal (U p (0) = 0),i.e., in the middle 
between two neighboring planes. So, for < y < 1, we have U p {y) = —U e (l — y) with U e (y) 
defined in flA.l|) . Note that only three of seven fitting parameters a^, yi,b, and c in ()A.1|) are 
independent. Using continuity conditions of the potential and corresponding electric field in 
points 2/1,2/2, we can, for example, express the rest parameters via y±, a\, and a 2 : 



a 2 - yjai(ai + a 2 ) y l (l - yx){a x + a 2 ) - 1 



°3 = z -7- 7^7 ; r, 2/2 



a x + (1 - 2/i) 2 (ai + a 2 ) ' a 2 - y 1 (o 1 + a 2 ) 



-2yi(ai + a 2 ) , c = 1 + ^(ai + a 2 ) 



(A.2) 



Independent fitting parameters along with £7 p z and d p j are listed in Table lA~T1 for some crystal 
planes. 



Table A.l: Potential well depths, U p i(eV), inter- planar distances, d p i(A), and parameters of 
the fit (|A.1|) for some crystal planes. 



Crystal, plane 


u pl 


d p i 


Vi 


CLi 


a 2 


diamond (110) 


23.54 


1.261 


0.109 


10.57 


1.994 


diamond (001) 


12.06 


0.892 


0.152 


7.16 


1.623 


silicon (110) 


21.27 


1.920 


0.100 


13.45 


2.750 


silicon (001) 


11.73 


1.358 


0.140 


8.63 


2.018 


copper (110) 


34.14 


1.278 


0.161 


7.00 


1.660 


iron (110) 


68.88 


2.027 


0.084 


17.38 


3.104 


tungsten (110) 


132.69 


2.238 


0.054 


34.90 


5.656 



The equation of one - dimensional motion in the potential U(y) is 

y = ±K^z-U(y)/U p i, k = 29 p i/ d pl , Q vi = ^2U p i/e 



(A.3) 



where z = s±/U p i = U(y)/U p i+ (y/n) 2 is the integral of motion (transverse energy in units of 
Upi). Remember now that the integration over time in formulas describing a radiation ( see 
Section 1) is performed at fixed value of z when, by means of a time shift, any initial conditions 
may be reduced to the standard one x(0) = providing the condition v x (— t) = v x {t). Then 
we have v n = v_ n = v* in the velocity fourier series 

00 

v x{t) = ^2 v n exp(inu t) , 



where c^o — 27r/T is the frequency of motion and T being its period. Finally, 

T/2 

v n = — dtv x (t) cos(nuj t) . 



(A.4) 



There is an additional symmetry, v x (t±T/2) = —v x (t), for channelled particles when Eq. ()A.4jl 
passes into 

l (~y)~T J dtv '^ t ) sin ( nuJ o t ) ■ ( A - 5 ) 





vf 1 = sin I 
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From this equation, even harmonics vanish at channelling (at —1 < z < for electrons 
and at < z < 1 for positrons). Let us define a quantity, v n , which is related to v n by 
v n = v n 8 p i/g(z) where g(z) = kT/4 at channelling and g(z) = k,T/2 at over-barrier motion. 
Solving Eq. (jA.3)l and taking elementary integrals in Eqs. (JA.4j) . (JA.5j) . we obtain for electrons 



9(z le - z) 



TC Z + 1 



4 V a>i 



1 



Ql + a 2 



Z — Z\ e 



cti + a 2 
Ql - ai 



\/z - z le Q e sin *i e 
(A.6) 



a%yi cos * le 



+ 9(z - z 2e 



a 3 — a 2 
Ql + 0,3 



y/z - Z 2e Qe Sin ^ 2e — a 3 (l — y 2 ) COS ^ 2 e 



where z le = U e {y x )/U pl = a x y\ - 1 , z 2e = U e (y 2 )/U pl = -a 3 (l - y 2 ) 2 , -1 < z u < z 2e < 



and 



7T 



C/ e (^) = 0(z le - z)— — + 9(z - Z\. 



0{z 2e - z) 



+ 6(z - z 2e 



arcsm 



l + zu 
l + z 



+ In 



am + y/ci2(z - zu) \ 
^\a 2 (z + c)-byA\) 



1 / Qij/i + A/a 2 (^ ~ Zle) 



«2 \a 3 (l - y 2 ) + y/a^z - z 2e ) 



+ 



(A.7) 



+ 



In 



(1 - y2)y/a^+ y/z - Z 2t 



Mr ^ 
Wie = — — arcsm 



1 + Zie 
l + Z 



2c 



Qc 



9e{Z) 



:lll 



(1 - y2)Jai+ y/z - Z 2c 



Qc 



29e{z, 



■[i + e(z)\ 



nuo, 



(e) 



ft 



Analogously, we find for positrons 



0« = ^-z) 



7T / £ 



<W + 5n,-l ) + 



4 V «3 

a 3 (l - 1/2) cos V lp + 0(z - z 2p ) 



Ql - a 2 

ai + a 2 
Ql + at 



9(z - z lp ) 



a 3 — a 2 



Ql - 0,3 



sfz - z lp Q p sin^ip 

(a.s; 

\Jz - z 2p Q p sin - a x yi cos ^ 2p 



where z\ p = —z 2e , z 2p = —z\ e , < z\ p < z 2p < 1 , and 



7T 



9{z 2p - z) 




Zip 



arcsm ( \ / — | + arccos 



\yJa 2 (z-c)+V/A 



+ 9{z - z 2p ) 



: arctan 



/ a\y\\Ja 2 {z - z lp ) - q 3 (l - y 2 )\J a 2 [z - z 2p ) 
V 01032/1(1 - 2/2) + a2A/( z ~ ^ipX^ - ^2 P ) 



+ 



(A.9) 



. 1 Zip 

—— arcsm * / — 

\/03 VV 2! 



2p 



2) . in' 



Tin 



2gJz i 



.[l + 9(z-l)) 



mo. 



(p) 
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As expected, differences in electron and positron motion diminish when z (transverse energy) 
increases. We obtain from above formulas for z ^> 1 (high above the potential barrier when 
z ~ (ip/Opi) 2 where ij) is the angle of a particle velocity w.r.t. the plane ) g e ~ g p ~ 
and therefore Uq ~ cjg — KKy/z — 27r0/dpj = qip. From Eq. ()A.6|l we find for z ^> 1 

1 

1 If 

Vn — 7 — ^[( a i + a 2) sin(7rra/i) + (a 3 - a 2 ) sin(7rm/ 2 )] = -Trrr / dy cos(7my)U e (y) . 
{7m) 6 2zUpi J 

o 

So, the velocity Fourier transform is expressed via Fourier coefficients of a potential. This 
happens when the rectilinear trajectory approximation is valid. If we substitute the original 
(periodic) plane potential for U e (y) in this integral, c4 takes the form dn^ ~ G(nq)/(2zU p i). 
Then the quantity jv n appearing in (|TU|) reads: jv n ~ G{nq) / {mip). At z 3> 1 , the velocity 
Fourier transform for positrons differs from that for electrons by the factor, — cos(7rn), which 
does not lead to some changes in (fTUj) as only v\ enters into this formula. 
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